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ABSTRACT 


attenuation  by  air  of  microwave  and  submillimeter  radiation. 
Besides  the  horizontal  attenuation#  the  vertical  attenuation  from 
various  levels  down  to  the  ground  and  out  into  space  is  calcula- 
ted for  a fixed  frequency.  The  line  profile  and  atmospheric  mo- 
del can  be  selected  from  among  several.  Comparison  is  made  with 
other  calculations,  and  with  experiments.  Possibilities  for 
improving  the  program  are  discussed. 
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TECHNICAL  REPORT  T-2/306-3-14 
PRESENT  STATUS  OF  THE  RRI  SLANT  - PATH 
ABSORPTION  MODEL  (SLAM)  COMPUTER  PROGRAM 

I.  INTRODUCTION 

A3  part  of  an  atmospheric  propagation  study  sponsored  by 
DARPA,  RRI  has  developed  a computer  program  to  calculate 
atmospheric  attenuation  in  the  microwave  and  submillimeter 
regions  of  the  electromagnetic  spectrum.  This  program  is 
intended  for  analyzing  communications  systems  and  other  appli- 
cations in  this  spectral  region  not  only  at  ground  level  but 
<also  with  transmitters  and  receivers  at  higher  altitudes. 

Hence,  emphasis  has  been  placed  on  calculating  the  total 
attenuation  down  to  the  ground  and  out  into  3pace  at  any 
given  frequency,  at  a set  of  reasonably  spaced  atmospheric 
levels.  The  program  is  written  in  FORTRAN1  and  may  be  looked 

upon  as  a modification  of  a program  originated  by  McClatchey, 
et.  al.  of  AFCRL. 2 

The  inputs  to  the  program  consist  of  an  atmospheric  model, 
a spectral  line  compilation,  and  a control  file  giving  the 
choice  of  parameters  to  be  used  in  the  calculations.  Currently 
six  atmospheric  models  are  available  for  use.  These  are  the 
"Tropical,"  "Midlatitude  Summer,"  "Midlatitude  Winter,"  "f ubartic 
Summer,"  "Subarctic  Winter,"  and  "U.S.  Standard  Atmosphere,  1962" 
all  as  formulated  by  McClatchey,  et.  al.3.  The  "Midlatitude  Winter 
model  is  reproduced  in  Table  I.  This  model  gives  the  pressure, 
temperature,  and  water  vapor  and  ozone  contents  as  a function  of 
altitude.  The  other  species  currently  represented  in  the  spectral 
line  compilation,  namely  oxygen  and  carbon  monoxide,  are  assumed 
to  have  constant  mixing  ratios.  Other  atmospheric  models  can  be 
used  by  changing  the  input  file.  The  spectral  line  compilation 
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TABLE  I 

MIDLATITUDE  WINTER  MODEL  USED  IN  SLAM  PROGRAM 


HEIGHT 

PRESSURE 

TEMP 

WATER 

OZONE 

(KM) 

(MBAR) 

«K) 

(Q/M#*3) 

(G/M**3) 

O* 

1 

* 0180E+03 

272.2 

3 

5 

E^OO 

6 

K> 

O 

■ 

U 1 

o 

o 

o 

o 

• 

1 • 

8 

973 

E*02 

268.7 

2 

5 

E*00 

6 

6 

E-05 

2* 

7 

897 

£♦02 

265.2 

1 

8 

E*00 

6 

9 

E-05 

3* 

6 

938 

E602 

241.7 

1 

2 

E^OO 

6 

9 

E-05 

6. 

6 

081 

E^02 

255.7 

6 

6 

E-01 

6 

9 

E-05 

5. 

5 

313 

E»02 

269.7 

3 

8 

E-01 

5 

8 

E-05 

4* 

6 

627 

E»02 

263.7 

2 

1 

E-01 

4 

6 

E-05 

7. 

6 

016 

E*02 

237.7 

8 

6 

E-02 

7 

7 

E-06 

8* 

3 

673 

E+02 

231.7 

3 

5 

E-02 

9 

0 

E-08 

9* 

2 

992 

E»02 

225.7 

1 

6 

E-02 

1 

2 

E-06 

10* 

2 

568 

E+02 

219.7 

7 

5 

E-03 

1 

6 

E-06 

11* 

2 

199 

E*02 

219.2 

6 

9 

E-03 

2 

1 

E-06 

12* 

882 

E»02 

218.7 

6 

0 

E-03 

2 

4 

E-06 

13. 

610 

E*02 

218.2 

1 

8 

E-03 

3 

0 

E-06 

16. 

378 

£♦02 

217.7 

1 

0 

E-03 

3 

2 

E-06 

IS. 

178 

E»02 

217.2 

7 

6 

E-06 

3 

6 

E-06 

16. 

007 

£♦02 

214.7 

6 

6 

E-06 

3 

6 

E-06 

17. 

610 

E*Oi 

216.2 

5 

6 

E-06 

3 

9 

E-06 

18. 

350 

£♦01 

215.7 

5 

0 

E-06 

6 

1 

E-06 

19. 

280 

E*Oi 

215.2 

6 

9 

E-06 

6 

3 

E-06 

20. 

5 

370 

E+Oi 

215.2 

6 

5 

E-06 

6 

5 

E-06 

21. 

580 

£♦01 

215.2 

5 

1 

E-06 

6 

3 

E-06 

22. 

910 

£♦01 

215.2 

5 

1 

E-06 

6 

3 

E-06 

23. 

360 

E^Ol 

216.2 

6 

6 

E-06 

3 

9 

E-06 

26. 

860 

E*Oi 

215.2 

6 

0 

E-06 

3 

6 

E-06 

25. 

630 

E + 01 

216.2 

6 

7 

E-06 

3 

6 

E-06 

30. 

110 

E»01 

217.6 

3 

4 

E-06 

1 

9 

E-06 

38. 

180 

£♦00 

227.8 

1 

1 

E-06 

9 

1 

E-06 

60. 

2 

530 

£♦00 

263.2 

6 

3 

E-05 

6 

1 

E-06 

65. 

290 

E+00 

258.5 

1 

9 

E-05 

1 

3 

E-06 

50. 

820 

E-01 

265.7 

6 

3 

E-06 

6 

3 

E-06 

70. 

6 

670 

E-02 

230.7 

1 

6 

E-07 

8 

6 

E-08 

100. 

3 

000 

E-06 

210.2 

1 

0 

E-09 

6 

3 

E — 1 1 

T-2/306-3-14 


RIVERSIDE  RESEARCH  INSTITUTE 


used  has  been  taken  in  large  part  from  the  first  section  of  a 
magnetic  tape  obtained  from  AFCRL4.  In  most  of  the  0 to  550  cm'"1 
spectral  region  under  consideration,  i.e.,  in  the  "pure  rotation" 
region  below  430  cm  , this  tape  includes  only  water,  ozone,  and  the 
microwave  spectrum  of  oxygen.  Consequently,  this  compilation 
has  been  modified  by  adding  the  results  recently  obtained  at 

C 

RRI  on  the  submillimeter  and  microwave  spectrum  of  oxygen  and 

on  the  rotational  spectrum  of  carbon  monoxide.8  (The  oxygen 

lines  were  added  for  three  isotopic  species:  18ol80, 180180  and 
18  18 

0 0.  Both  the  ground  vibrational  level  and  the  first  vibra- 

tionally  excited  state  were  taken  into  accound  for  l60160.) 

A sample  of  the  compilation  is  shown  in  Table  II.  The  first 
column  gives  the  frequency  in  cm-1?  the  second  column,  the 
(modified)  integrated  line  strength  at  296K,  in  cm-1/molecules  cm-2 
(modification  described  below) ; the  third  column,  the  pressure 
coefficient  of  the  line  half-width  in  cm-1/atm;  the  fourth  column, 
the  energy  of  the  lower  state  in  cm-1;  the  fifth  column,  the  date 
(month  and  year)  of  insertion  into  the  file;  the  sixth  column, 
the  isotope  concerned  (e.e.,  68  is  0l6018) ; and  the  seventh  column, 
the  molecular  constituent,  the  integers  M=l,3,5,7  standing  for 
water  vapor,  ozone,  carbon  monoxide,  and  oxygen,  respectively. 

The  data  on  the  tape  received  from  AFCRL  was  first  con- 
verted to  the  9-track  format  required  by  the  Xerox  Sigma  9 com- 
puter and  then  was  placed  into  a magnetic  disc  file  and  converted 
from  BCD  to  binary  format  (for  increased  speed) . The  control 
file  (Table  ill)  contains  information  on  the  first  and  ladt  fre- 
quency to  be  used,  the  frequency  increment,  the  quantity  BOUND 
(which  gives  the  frequency  range  within  which  lines  are  to  be 
summed  to  obtain  the  attenuation  at  a fixed  frequency) , the 
choice  of  the  spectral  line  profile  to  be  used,  and  a choice  of: 

(a)  summing  over  all  the  molecular  constituents  to  get  the  total 
attenuation  (zero  in  the  sixth  column) , (b)  just  considering  the 
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TABLE  II 

Sample  of  Spectral  Line  Compilation 


FRfcO 

strength 

half 

wIOTm 

NO  OF  RElUHDS- 

40 

•3u7 

•121E-22 

• 1100 

• 3**  1 

• 1996-22 

• 1100 

•3c9 

•199E-22 

• 1100 

• 496 

• 1836-22 

• 1100 

• 539 

•1576-22 

• 1100 

• 742 

•121E-P1 

• 0811 

• 7*6 

•355E-22 

• 1100 

• 856 

•192E-22 

• 1100 

• 966 

•3266-22 

• 1100 

1 • 0 j2 

•6296-22 

• 1100 

1 • Uy7 

.*896-22 

• 1100 

1-018 

•169E-22 

• 1100 

1 • 2o2 

•395E-22 

• 1100 

1 • ?o2 

•683E-22 

.1100 

1*4*9 

•254E-22 

• 1100 

1 • 4*6 

•492E-22 

*1100 

1 -463 

• UOE-22 

*1100 

1*4*7 

•1756-22 

• 1100 

1*650 

*5786-2/ 

• 0320 

1 »6t»7 

•1716-26 

.0320 

l*Ce9 

•3186-22 

• 1100 

1*6/6 

•6796-24 

• 0900 

1*664 

•476E-26 

.0320 

1 • 7jl 

•125E-25 

.0320 

1 • 7 j3 

•2126-22 

• 1100 

1*718 

•3106-25 

.0320 

1*7J4 

•575E-22 

• 1100 

1*735 

• 725E-25 

.0320 

1*753 

•1606-26 

• 0320 

1*76** 

•437E-27 

• 0320 

1*7/0 

•3326-26 

• 0320 

1*7/3 

•621E-27 

• 0320 

1 *7*1 

•868E-27 

• 0320 

1 *7«8 

*6506-26 

• 0320 

1*769 

*1206-26 

• 0320 

1*7*1 

•7B6E-22 

• 1100 

1*797 

• 1636-26 

• 0320 

1 *8U0 

*6356-27 

• 0380 

1*806 

•1206-23 

• 0380 

1 • Bu6 

•2186-26 

• 0350 

NO  OF  RECORDS- 

40 

1*816 

•2896-26 

.0380 

1*819 

• 1106-26 

• 0350 

1*8*3 

•3766-26 

• 0370 

1*8*6 

•207E-23 

• 0350 

1*831 

•4836-26 

• 0350 

1*839 

* 1 77E-26 

• 0370 

1*840 

•6106-26 

.0360 

1*842 

•3366-23 

.0370 

1*846 

•4816-22 

• 1100 

ENEHGY  DATE  ISO  CONSTITUENT 
TOPE 


205*328 

102  < 

666 

3 

50*302 

102 

666 

3 

8*022 

102 

666 

3 

281*833 

102 

666 

3 

365.692 

102 

666 

3 

66*. 512 

23 

161 

1 

171*501 

102 

666 

3 

127*266 

102 

666 

3 

3014671 

102 

666 

3 

128.119 

102 

666 

3 

100*572 

102 

666 

3 

157*167 

102 

666 

3 

261.831 

102 

666 

3 

156*903 

102 

666 

3 

2*519 

102 

666 

3 

77*082 

102 

666 

3 

702*316 

102 

666 

3 

190*212 

102 

666 

3 

2460*776 

75 

66 

7 

2230*625 

75 

66 

7 

667*787 

102 

666 

3 

155*389 

52 

162 

1 

2011.215 

76 

66 

7 

1803.180 

75 

66 

7 

579*061 

102 

666 

3 

1606*353 

75 

66 

7 

323*620 

102 

666 

3 

1420*767 

75 

66 

7 

1246*652 

75 

66 

7 

1178*121 

75 

68 

7 

1083*636 

75 

66 

7 

1099*777 

75 

68 

7 

1026*107 

75 

68 

7 

931*765 

75 

66 

7 

951*113 

75 

68 

7 

34*251 

102 

666 

3 

•80*799 

75 

68 

7 

2339*133 

75 

66 

7 

791*605 

75 

66 

7 

813*167 

75 

68 

7 

761*219 

75 

61 

7 

2211*583 

75 

66 

7 

685*959 

75 

68 

7 

662*637 

75 

66 

7 

626*3l> 

75 

68 

7 

2095*301 

75 

i 66 

7 

569*509 

75 

68 

7 

564 i 863 

75 

66 

7 

281*833 

102 

1 666 

3 
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attenuation  due  to  one  particular  molecular  species  (M  in  the 
sixth  column),  or  (c)  considering  the  attenuation  due  to  all 
species  except  one  (-M  in  the  sixth  column) . The  number  in  the 
seventh  column  chooses  the  caption  to  go  with  the  atmospheric 
models  described  above.  Three  choices  of  line  profile  are 
currently  available t Lorentz,  van  Vleck-Weisskopf  and  Gross 
(kinetic).  ' ' These  are  denoted  by  1,2,  or  3 respectively  in 
the  fifth  column.  When  the  frequency  at  which  the  attenuation 
is  being  calculated  coincides  within  .000S  cm-1  with  one  of  the 
spectral  line  frequencies  in  the  data  file,  the  three  profiles 
mentioned  have  been  replaced  by  a Voigt  profile  for  the  parti- 
cular line  for  which  this  is  true,  in  this  case  the  center  of 
the  line  is  assumed  to  coincide  exactly  with  the  frequency  at 
which  the  attenuation  is  being  calculated  and  the  Voigt  profile 
at  line  center  is  used.  This  was  found  necessary  because  the 
vertical  attenuation  from  a given  high  altitude  out  into  space 
is  greatly  overestimated  unless  line  broadening  due  to  the 
Doppler  effect  is  included.  The  Voigt  profile  includes  a combi- 
nation of  collisional  and  Doppler  broadening.  The  choice  of 
Voigt  profile  at  line  center  was  made  because  this  gives  the 
maximum  attenuation  of  which  the  line  is  capable.  Any  other 
frequency  off  line  center  will  give  a lower  attenuation.  The 
SLAM  program  currently  in  use  is  shown  in  the  Appendix. 
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TABLE  III 


EXAMPLE  OF  CONTROL  FILE 


INITIAL 

FREQ* 

FINAL 

FREQ* 

FREQ* 

INTERVAL 

BOUND 

LINE 

shape 

CON* 

ATM* 

MODEL 

2*0*0 

2*040 

'001 

5* 

1 

7 

6 

-6- 
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II . METHODS  OF  CALCULATION 


The  atmospheric  model  directly  gives  the  concentrations 
m grams  pec  cubic  meter  for  water  vapor.  Cl (N) , and  for 
ozone , C3U),  at  the  altitude  represented  by  the  (sequential) 
integer,  N.  These  mass  densities  are  converted  into  the 
number  of  molecules  in  a column  1 km  long  by  1 cm2  cross- 
sectional  area,  denoted  by  W(1,N)  and  W(3,N),  repectively, 
by  the  formulae : * 

W(1,N)  = Cl  (N)  * 3 . 346E+21 
W (3 ,N)  = C3(N)  * 1.2546E+21. 

For  those  molecules,  M,  which  have  a constant  mixing 

ratio,  AM(M),  the  quantities  W(M,N)  are  computed  by  the 
formula : 

W(M,N)  = ( . 724270E+24) *P (N) *AM (M)/T (N) . 

where  the  numerical  factor  is  the  reciprocal  of  k,  the 
Boltzmann  constant, times  a power  of  10.  P(N)  and  T(N)  are 
respectively,  the  pressure  in  millibars  and  the  temperature 
in  degrees  Kelvin,  as  given  by  the  atmospheric  model.  M = 5 
and  7 correspond  to  carbon  monoxide  and  oxygen,  respectively. 
The  mixing  ratio  AM(M)  represents  the  part  of  the  atmosphere 

by  Volume  that  consists  of  the  naturally  occurring  isotopic 
mixture  of  constituent  M: 

AM(5)  = . 075E-06 
AM(7)  = .2095 


* W(M,N) 
M(12c)  = 


= CM (N)/ (10* (Molecular  Weight) * (1/12*M  12C 
mass  of  carbon-twelve  atom. 


)), 
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Next,  the  population  of  each  quantum  state  capable  of  absorbing 
energy  in  the  spectral  region  under  consideration  must  be 
determined.  The  fractional  population  is  given  by  the  expression 

EXP  (-eAt)/Q{t)  , 

the  Boltzmann  exponential  population  factor  EXP(-E/kT)  divided 
by  the  partition  function  Q (T) , at  the  temperature  T.  In 
absorption,  E is  the  energy  of  the  lower  state;  k is  Boltzmann's 
constant.  The  exponential  population  factor  and  the  partition 
function  have  been  evaluated  at  the  standard  temperature  T0=296K 
and  the  results  are  included  in  the  spectral  line  data  file  as  a 
factor  of  the  line  strength.  The  corrections  to  these  factors 

necessary  at  temperatures  T other  than  TO  are  computed  by  the  program. 

I 4-6 

The  exact  quantum  partition  function  was  evaluated  at  TC«*296K,  and 

the  temperature  correction  is  approximated  by  assumimg  the  temperature 

dependency  to  be  that  of  the  classical  rotational  partition 

function  (see  e.g..  Appendix  D of  Ref.  5).  Thus,  linear 

molecules  are  assumed  to  have  a rotational  function  that  varies 

directly  as  the  temperature,  and  non-linear  molecules  one 

chat  varies  as  the  3/2 -power  of  the  temperature.  The 

vibrational  partition  function  at  the  standard  temperature 

is  included  but  its  temperature  dependence  is  neglected.  The 

quantity  necessary  to  adjust  the  exponential  population  factor 

from  temperature  TO  to  T (N)  is 

CS1(N)  = (TO-T (N) ) / (TO*T (N) * . 6946) 

and  the  quantity  necessary  to  adjust  the  partition  function 
from  TO  to  T (N)  is  given  by 

CS2(M,N)  = TO/T (N) 
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for  linear  molecules,  and  by 

CS2(M,N)  = (TO/T (N) ) **  1.5 
for  non-linear  molecules. 

The  stimulated  emission  at  TO  was  included  in  the  line 
strength  on  the  original  AFCRL  data  tape.  However,  the  LBL 
program  listed  in  Ref.  (2)  did  not  include  corrections  to  the 
stimulated-emission  factor  at  other  temperatures.  In  the 
present  SLAM  program  this  stimulated-emission  factor  has  been 
divided  out  of  the  line  strengths  in  the  data  file  and  the 
SLAM  program  itself  computes  the  required  stimulated-emission 
factor  at  the  arbitrary  temperature  T(N).  This  factor  is 

1.  -EXP (-4 . 86378E-03*GNU (I) ) , 

where  GNU (I)  is  the  wavenumber  (cm-1)  of  the  spectral  line 
whose  contribution  to  the  attenuation  we  are  considering.  The 
final  value  of  the  line  strength,  when  multiplied  by  the 
exponential  population  factor  and  the  stimulated  emission 
factor,  and  then  divided  by  the  partition  function  is 
there fore : 

SS=S (I ) *CS2 (M,N)*EXP (-EPD (l)*CSl (N) ) 

* (1 . -EXP (-4 . 86378E-03*GNU (I) ) ) . 

Here  S(I)  is  the  line  strength  as  it  appears  in  the  data  file 
(DATA  2)  and  EPD(I)  is  the  energy  of  the  lower  state  from  which 
the  absorption  arises.  S(I)  also  includes  as  a factor  the 
isotopic  abundance  ratios  of  molecules  such  as  Cf^O1®  which 
are  included  in  the  file. 

Next,  the  line  width  must  be  calculated.  The  data  file 
contains  the  half -width  of  the  line  at  half -maximum  at  the 
standard  temperature  TO,  and  pressure  PO  = 1013  millibars. 

When  the  linewidth  is  assumed  to  be  proportional  to  the 
pressure,  the  temperature  and  pressure  dependence  of  the  line 
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width  may  be  taken  into  account  by  calculating  the  factor: 
CA(N)  = ((T0/T(N)  )**0.5)*(P(N)/P0)  . 

This  factor,  which  is  a first  approximation,  assumes  that  a 1 
lines  of  all  species  have  the  same  temperature  dependence  as 
would  occur  when  the  collision  diameters  are  independent  of 
temperature.  The  half-width  is  therefore  given  by 

ALPHA1  = ALPHA(I)*CA(N)  . 

where  ALPHA (I)  is  the  half  width  at  P0,T0  of  the  Ith  line 
given  in  the  file. 

Three  line  shapes  can  be  chosen  by  the  user  for  the 
calculations:  the  Lorentz,  van  Vleck-Weisskopf  (wwf),  and 
"kinetic"  (Gross)  line  shapes.  For  any  line  shape,  the 
auxiliary  quantitities 

Z = V-GNU(I) 

Zl'-'V+GNU  (I ) 

are  calculated  first,  where  V is  the  wavenumber  at  which  the 

attenuation  is  to  be  calculated,  and  GNU (I)  is  the  wavenumber 
of  the  Ith  line.  If 

ABS (Z) *LE* .0005 

the  Voigt  profile10  is  used  for  the  calculation. 

For  the  Lorentz  line  shape,  the  quantity 

SUM1 (M) =SS*ALPHA1/ (Z**2+ALPHA1**2 ) 

is  calculated. 

The  Van  Vleck-Weisskopf  shape  requires  first  the  cal- 
culation of 

WWF=(1‘/(Z**2+ALPHA1**2) ) + (1  */(Zl**2+ALPHAl**2 ) ) 
and  then  the  calculation  of 

SUM1  (M)  = (SS*ALPHA1*  (V**2 )/ (GNU  (I )**  2 ) )*WWF. 


-10- 


T-2/306-3-14 


nlftnolut  HtitAhlih  INSIIIUIt 


The  Gross  (kinetic)  line  shape  is  obtained  by  first 
calculating  the  factor 

GR=(  (Z*Z1)**2)  + (4.*  (GNU  (l)**2)*  (ALPHA1**2) ) 
and  then  calculating 

SUM1  (M)=SS*ALPHA1*4«*(GNU(I)**2)/GR. 

The  Doppler  broadening  is  determined  by  the  parameter  BETA 
which  is  the  half-width  divided  by  J In  2 of  the  Doppler  broadened 
line  at  half-maximum.  The  quantities  DOP(M)  which  depend  on 
the  masses  of  the  individual  molecular  species  M are  stored 
in  the  program  as  DATA  and  are  used  to  compute : 

D0P1 (Nl ,N2 ) =DOP (Nl) *T (N2 ) ** . 5 

which  now  depends  on  the  altitude  and  atmospheric  model  through 
dependence  on  the  temperature  T(N2).  Then  BETA  is  given  by: 

BETA=GNU (I ) *D0P1 (M, N) 

where  the  dependence  on  the  line  frequency  is  now  included. 

The  shape  of  the  Voigt  profile  depends  on  the  dimensionless 
variable  Y: 

Y=ALPHA 1/BETA , 

which  is  the  ratio  of  the  Lorentz  half-width  to  the  Doppler 
parameter. 

In  the  special  case  of  the  center  of  a line,  calculation 
of  the  Voigt  profile  can  be  reduced  to  the  calculation  of  the 
probability  integral.  When  GNU(I)=30  cm-1  and  Z=.001  cm-1 
the  parameters  are  already  such  that  the  Voigt  profile  can  be 
approximated  by  a Lorentz  shape.  Therefore  in  the  first  approxi- 
mation it  was  not  thought  necessary  to  replace  the  other  profiles 
by  a Voigt  one  except  possibly  for  the  line  nearest  the  frequency 
at  which  the  attenuation  is  being  calculated.  The  expression  to 
be  evaluated  is: 

2 oY2f  -q2 

BETA  S »'  e dq 
Y 

where  a factor  1/tt  has  been  omitted. 
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When  (Y.GE.5.)  the  asymptotic  formula  for  the  probability  inte- 
gral is  used, giving: 

SUM1 (M) =SS* ( 1 . /Y- . 5/Y3+ . 75/Y5 

-1 . 875/Y7+6 . 5625/Y9) /BETA 

where  Y3  is  the  third  power  of  Y,  etc.  When  Y is  less  than  5 an 
approximation  of  Hastings^  is  used.  Let: 

TY=1./(1.+.3275911*Y) 

and  TY3  be  the  third  power  of  TY,  etc.  Then  use  of  Hastings' 
approximation  gives : 

SUMl (M) =SS*1. 772454* ( . 2548296*TY 

- .2844967*TY2+1.421414*TY3-1.453152*TY4 

+ 1 . 061405*TY5 ) /BETA . 

All  line  shapes  given  above  neglect  a possible  shift  of 

resonant  frequency  with  pressure  of  about  .01  cm-1  per  atmos- 
12 

phere . This  sets  a limit  on  the  precision  with  which  it  is 
desirable  to  give  the  resonant  wavenumbers  in  the  data  file. 

(Note  that  .001  cm  1=30  MHz,  a reasonable  modulation  bandwidth.) 
In  addition,  the  exact  shape  of  the  spectral  lines  far  from  the 
resonant  frequencies  is  not  settled,  especially  for  the  case  of 
water  vapor.  In  this  case,  besides  the  usual  monomer  form, 
water  dimers  may  exist  and  contribute  appreciably  to  the  absorp- 
tion. 

Currently,  the  continuum  contributions  (due  primarily  to  Nj 
and  in  part  to  non-re sonant  absorption  by  Oj)  have  not  yet  been 
included  in  the  SLAM  program.  Thus,  the  attenuation  in  the 
vicinity  of  absorption  minima  will  be  somewhat  underestimated. 

To  determj  '.e  the  attenuation,  the  quantities  SUMl(M)  from 
the  different  lines  of  the  species  M must  be  added  at  a fixed 
wavenumber : 

CAYl  (M)  =CAY1  (M I +SUM1  (M)  , 

where  this  sum  is  iterated  over  all  the  relevant  lines. 
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T^e  quantity  CAYl(M)  is  next  multiplied  by  W(M,N)  and  the 
result  summed  over  M: 

CAY=CAY+CAYl(M)*W(M,N) . 

Finally,  a normalization  factor  of  1/if  omitted  in  the 
definition  of  the  lineshape,  and  conversion  of  units  into 
dBAm,  gives  the  formula 

OPD  (IV)  =CAY*  1 . 38246 

where  OPD  is  the  horizontal  attenuation  in  dBAm. 

In  addition  to  the  horizontal  attenuations  at  the  levels 
represented  by  N,  the  integrated  attenuations  from  these 
levels  down  to  the  ground  and  up  into  space  are  calculated. 
The  integration  is  performed  by  means  of  the  trapezoidal  rule 
An  iteration  of  the  formula 

SUM=SUM+. 5 * (H (N) -H (N-l) )* (OPD (IV) +SAVE) 

gives  the  attenuation  from  the  height  H(N)  down  to  the  ground 
level,  where  OPD (IV)  is  the  horizontal  attenuation  at  H(N), 
SAVE  is  the  horizontal  attenuation  at  the  height  H(N-l)  and 
SUM  is  initialized  to  zero  at  the  ground  level.  Then 

DOWN (N, IV) =SUM 

is  the  attenuation  in  dB  from  the  level  H(N)  down  to  the 
ground  . 

HOR  (N,  IV)  =0 PD  (IV) 

is  the  horizontal  attenuation  in  dBAm,  while 
UP  (N,  IV)  =SUM-D0WN  (N,  IV) 

iQ  attQnuation  in  dB  from  the  height  H (N)  out  into  space. 
In  the  latter  formula  SUM  represents  the  attenuation  from 
the  highest  level  down  to  the  ground.  The  wavenumber  is 
then  incremented  and  the  calculation  begun  anew. 
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Typical  output  in  tabular  form  is  shown  in  Tables  IV  and  V. 
Columns  labeled  HEIGHT  give  the  altitude  in  kilometers  above 
mean  sea  level,  those  marked  HOR  give  the  attenuation  in  the 
horizontal  direction  in  dBAm,  those  marked  DOWN  give  the  total 
attenuation  in  dB  from  the  height  indicated  down  to  sea  level, 
while  those  marked  UP  give  the  total  attenuation  in  dB  from  the 
height  indicated  vertically  outwards  into  space.  The  numbers 
marked  FREQUENCY  give  the  frequency  (in  cm"1)  at  which  the  cal- 
culations were  done.  From  the  heading,  we  note  that  a Van  Vleck- 
Weisskopf  line  profile  was  used  in  Table  IV,  that  BOUND  was 
20  cm  , that  all  constituents  on  the  data  file  were  used  in  the 
calculation  (shown  by  the  value  0 after  CONSTITUENT),  and  that 
the  U.S.  Standard  Atmosphere,  1962  was  used.  This  frequency  is 
in  a valley  between  two  water  vapor  absorption  lines.  Table  V 
shows  the  results  of  a calculation  identical  to  that  of  Table  iv, 
except  that  a Gross  profile  was  used  instead  of  the  Van  Vleck- 
Weisskopf  profile.  The  two  tesults  for  HOR  differ  by  about  25% 
at  sea  level,  but  the  difference  grows  considerably  larger  at 
high  altitudes.  At  the  peak  of  absorption  lines  little  difference 
is  found  in  the  calculated  attenuation  as  a function  of  line 
shape . 

Figures  1 through  7 show  typical  graphical  output.  Figures 
1 and  2 present  graphs  of  HOR,  DOWN  and  UP  corresponding  to  the 
numerical  output  of  Tables  IV  and  V.  The  caption  identifies  the 
various  curves  and  shows  whether  they  are  measured  in  dB  or 
dBAm,  i.e.,  whether  the  scale  of  ordinates  on  the  left  or  on 
the  right  is  to  be  used.  The  abscissa  is  the  HEIGHT,  and  the 
caption  gives  information  similar  to  that  in  the  tabular  output. 
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concerning  the  choice  of  parameters  for  the  calculation.  Com- 
parison of  Fig.  1 and  2 again  shows  the  effect  of  line  shape. 

Figures  3 and  4 plot  the  computed  attenuation  for  a fre- 
quency of  2.040  cm  \ which  is  at  the  peak  of  the  60  GHz  oxygen 
microwave  absorption  band.  Figure  3 shows  the  attenuation  com- 
puted as  described  above  with  the  Voigt  profile  substituted  for 
the  Lorentz  profile  for  the  nearest  line.  In  Figure  4 the 
Lorentz  profile  at  the  center  of  a line  has  been  used  for  this 
line  too.  Corresponding  numerical  values  are  shown  in  Tables  VI 
and  VII.  It  can  be  seen  that  HOR  begins  to  differ  appreciably 
for  the  two  cases  above  50  km.  This  lends  to  a considerable 
difference  in  the  computed  value  of  UP  at  all  altitudes.  Although 

the  Lorentz  line  shape  was  used  for  this  calculation  in  order  to 

13 

facilitate  comparison  with  the  work  of  Liebe  and  Welch,  the 
Van  Vleck-Weisskopf  profile  is  generally  considered  more  accurate 
in  the  microwave  region  and  the  Gross  profile  is  perferred  in  the 
submillimeter  region.  Comparison  of  Fig.  3 with  Fig.  1 or  2 
shows  considerable  difference  in  the  shape  of  the  altitude- 
dependent  curves.  The  relatively  fast  drop-off  of  the  HOR  curve 
in  Figs.  1 and  2 is  undoubtedly  due  to  the  small  scale-height 
for  water  vapor,  which  contributes  most  of  the  absorption  in 
the  troughs.  The  scale -height  for  oxygen,  however,  is  the 
same  as  that  of  the  atmosphere  as  a whole,  so  that  HOR  at 
the  peak  of  an  oxygen  absorption  line  does  not  fall  off  as 
rapidly. 

Figures  5,6,  and  7 are  calculated  at  a frequency  of 
14.169  cm  , which  is  the  peak  of  one  of  the  strong  molecular 
oxygen  submillimeter  lines.  Figure  5 is  calculated  with  all 
the  constituents  in  the  data  file  (CONSTITUENT  =0),  Fig.  6 with 
only  the  oxygen  lines  (CONSTITUENT  =7),  and  Fig.  7 with  all 
lines  except  those  of  oxygen  (CONSTITUENT  =-7)  . Otherwise, 
the  calculations  for  these  figures  were  done  using  identical 
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used  for  all  lines  except  the  center  line  for 
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fa lues  of  the  parameters.  Comparison  of  Figs.  5 and  6 shows 
that  the  general  altitude -dependence  of  the  attenuation  can  be 
traced  to  the  oxygen.  Figure  7 shows  that  the  other  consti- 
tuents make  an  appreciable  contribution  at  low  altitudes.  Thus 
the  importance  of  including  the  oxygen  submillimeter  spectrum 
in  this  spectral  region  is  illustrated. 
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IV.  COMPARISON  WITH  CALCULATIONS  BY  OTHERS 


Table  VIII  shows  the  results  of  varying  BOUND  within  the 
SLAM  program  for  a frequency  of  1.000  cm"1.  Also  shown  are  a 
comparison  with  calculations  by  Van  Vleck.14  in  these  calcula- 
tions only  lines  in  the  spectrum  of  water  vapor  are  taken  into 
account,  all  line  half -widths  are  put  equal  to  .1  cm"1  and  the 
temperature  and  pressure  shown  in  the  table  are  used  instead 
of  one  of  the  standard  atmospheric  models.  This  was  done  to 
make  a closer  comparison  with  Van  Vleck' s results.  The  quantity 
given  in  the  table  is  the  horizontal  attenuation  in  dBAm  per 
unit  density  of  water  vapor,  the  latter  expressed  in  g/m3. 

The  results  agree  to  about  3.6%  which  is  considered  satisfactory. 

Table  IX  shows  a comparison  of  the  results  of  calculations 
performed  by  SLAM  with  those  performed  at  AFCRL4.  The  results 
of  varying  BOUND  within  the  SLAM  program  for  this  higher  fre- 
quency are  also  shown.  The  discrepency  between  the  SLAM  and 
AFCRL  calculations  is  about  ten  percent.  This  is  considered  a 
reasonable  agreement;  further  comments  are  difficult  to  make 

since  full  documentation  concerning  the  AFCRL  calculations  are 
not  available  at  RRI. 

Table compares  the  results  of  SLAM  and  those  of  Liebe 
and  Welch  , using  the  Lorentz  profile  and  U.S.  Standard  Atmos- 
phere, 1962  used  by  those  investigators.  The  quantities  BOUND 
and  CONSTITUENT  are  set  so  that  only  the  oxygen  microwave  spec- 
trum contributes  to  the  attenuation.  The  discrepancy  between 
the  two  calculations  ranges  up  to  a factor  of  two  at  10  km.  it 
is  believed  that  this  discrepancy  can  be  traced  to  different 
assumptions  made  concerning  the  variation  of  line  half-width 
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FREQUENCY  - 29.7l3cnf'  VAN  VLECK-WEISSKOPF  PROFILE  MIDLATITUDE  WINTER  MODEL 

CONSTITUENT  * 0 
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with  altitude.  The  SLAM  program  assumes  line  half-widths  chara- 
cteristic of  well  separated,  isolated  lines,  which  is  justified 
for  low  pressures  of  the  broadening  gas.  However,  due  to  inter- 
action between  overlapping  lines,  at  higher  pressures  the  attenu- 
ation far  from  the  oxygen  band  is  overestimated  by  this  procedure. 
An  empirical  correction  is  often  made  by  decreasing  the  assumed 
line  half-widths  until  the  calculated  attenuation  outside  the 

band  matches  the  observed  result.  In  the  calculations  of  Liebe 
13 

and  Welch  this  empirically  corrected  line  half-width  has  been 
used  up  to  an  altitude  of  about  15  km  after  which  a gradual 
transition  to  the  line  half -widths  characteristic  of  the  isolated 
lines  is  made.  Since  at  10  km  the  overlapping  of  the  lines  has 
decreased  considerably  this  could  account  for  the  large  discre- 
pancy between  the  two  calculations  at  that  altitude. 
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V.  COMPARISON  WITH  EXPERIMENTS 


Table  XI  shows  a comparison  of  SLAM  output  with  experimental 
results  of  Becker  and  Autler15  and  of  Burroughs,  Jones  and 
Gebbie  . Only  attenuation  by  water  vapor  is  included  and  the 
results  are  normalized  to  1 g/m^  of  water  vapor.  Again,  instead 
of  one  of  the  standard  atmospheric  models,  the  temperature  and 
pressures  shown,  appropriate  for  a comparison  with  the  corres- 
ponding experiment,  have  been  used.  There  is  substantial  agree- 
ment between  experiment  and  the  results  of  either  the  Van  Vleck- 
Weisskopf  or  G-v-ss  profiles  at  .784  cm"1,  which  is  near  the 
22  GHz  water  vapor  absorption  line.  Further  from  this  line  the 
discrepancy  increases  to  a factor  of  about  50%  at  1.160  cm-1  and 
100%  at  1.340  cm  . The  Lorentz  profile  results  are  considerably 
more  discordant.  Tho  values  quoted  from  Burrou7hs,  Jones  and 
Gebbie  have  been  obtained  by  evaluating  their  least  squares 
fitted  polynomial,  which  is  a function  of  the  broadener  pressure, 
at  the  appropriate  value  of  the  atmospheric  pressure,  in  this 
evaluation  the  term  independent  of  the  broadener  pressure  has  been 
omitv.ed  as  being  due  to  the  broadening  of  the  water  vapor  lines 
by  water  vapor,  which  is  not  included  in  the  SLAM  program,  it 
can  be  seen  that  the  experimentally  determined  attenuations  are 
still  considerably  higher  than  the  calculated  ones. 

Table  XII  compares  the  SLAM  output  with  atmospheric  obser- 
vations by  Lo,  Fannin  and  Straiton17.  These  numbers  have  been 
taken  from  one  of  their  fitted  curves.  Also  shown  are  the  de- 
composition of  the  calculated  absorption  into  the  sum  of  absorp-’ 
tion  by  oxygen  and  by  water  vapor.  From  the  above  discussions 
it  follows  that  two  compensating  errors  are  present  in  these 
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TABLE  XTT 
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calculation*.  The  excess  attenuation  by  water  vapor  is  not  in- 
cluded in  SLAM,  while  the  attenuation  due  to  oxygen  is  overe5;tim- 
ated  due  to  neglect  of  the  interacting  line  mechanism.  Once 
again  the  Lorentz  profile  results  are  very  large.  Table  XIII 
shows  the  considerable  variation  of  the  Lorentz  profile  results 
with  BOUND.  These  results  can  be  attributed  to  the  large  low 
frequency  tail  from  the  Lorentz  profile  of  lines  at  high  fre- 
quencies . 
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VI.  FURTHER  REFINEMENTS  REQUIRED 


Several  items  which  were  not  yet  included  in  the  program 

and  which  might  have  an  appreciable  effect  on  the  attenuation 

are  discussed  below.  In  the  SLAM  program  so  far,  we  have 

neglected  the  continuum  attenuation  of  nitrogen.  Using  the 

18 

data  of  Gebbie,  et.  al.  we  estimate  the  horizontal  attenuation 
due  to  this  cause  at  ground  level  to  be  about  1.5  dB/km  at  a 
frequency  of  100  cm  the  peak  of  the  nitrogen  continuum  absorp- 
tion curve,  and  0.2  dB/km  at  a frequency  of  35  cm-1.  The  nitro- 
gen continuum  should  be  included,  but  it  is  not  likely  to  change 
the  general  shape  of  the  attenuation  curves.  Several  trace 
species  having  active  rotational  spectra  are  present  in  the 
stratosphere,  and  it  is  desirable  to  estimate  the  strengths  of 
their  absorption  lines.  These  include  nitrous  oxide,  nitric 
acid  vapor,  nitrogen  dioxide,  sulfur  dioxide,  and  nitric  oxide.19 

A low-altitude  effect  that  should  be  dealt  with  is  the 
excess  attenuation  of  water,  beyond  that  calculated  from  the 
individual  lines  of  the  water  monomer.  This  is  especially 
noticeable  in  the  troughs  between  the  peaks  of  the  water  monomer 
absorption.  This  excess  absorption  can  have  three  causes.  The 
line  profiles  of  the  water  monomer  may  be  different  from  that 
assumed,  leading  to  greater  absorption  in  the  wings  of  ihe  lines, 
and  hence  in  the  troughs.  Water  dimers29-23  can  contribute 
their  own  spectrum  of  absorption,  and  lastly,  short-lived  colli- 
sion complexes  of  the  water  molecule  with  the  broadening  molecule 
may  induce  a temporary  electric  dipole  moment,  leading  to  in- 
creased absorption.  These  three  possible  causes  can  be  largely 
separated  from  each  other  by  considering  how  they  scale  with  the 
physical  variables  of  the  problem:  water-vapor  pressure,  pressure 
of  the  external  broadening  gas,  and  temperature. 
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Various  cases  can  be  distinguished:  if  S is  the  integrated 
strength  of  a line  and  a its  half-width  at  ha If -maximum,  then  at 
the  center  of  an  isolated  line,  the  attenuation  is  proportional 
to  S/a;  in  the  wings,  it  is  proportional  tc  S a;  and  for  a group 
of  lines  with  strength  S that  overlap  due  to  the  collisional 
broadening,  it  is  S.  For  water  monomers,  S is  proportional  to 
the  partial  pressure  of  the  water  vapor;  for  dimers,  it  is  pro- 
portional to  the  square  of  the  partial  pressure,  and  for  colli- 
sion complexes  of  water  molecules  with  nitrogen,  it  is  propor'- 
tional  to  the  product  of  the  partial  pressures  of  water  vapor  and 
nitrogen.  Likewise,  the  quantity  a is  proportional  to  the  partial 
pressures  of  nitrogen  and  water  vapor  for  broadening  by  nitrogen 
and  water  vapor,  respectively.  For  collisionally  induced  absorp- 
tion, a is  independent  of  any  of  the  pressures.  From  the  data 
of  Burroughs,  Jones,  and  Gebbie,16  taken  at  311  and  337  Mm  in 
pure  water  vapor,  the  attenuation  is  proportional  to  the  square 
of  the  water  vapor  pressure.  From  the  above  scaling  laws,  this 
might  be  due  either  to  attenuation  from  the  wings  of  well-se- 
parated monomer  lines  or  to  overlapping  dimer  lines.  However, 
the  rapid  change  of  attenuation  with  temperature  was  interpreted 
by  these  experimenters  as  evidence  of  water  dimers. 

At  altitudes  higher  than  about  40  km,  Doppler  broadening 
will  be  important,  its  neglect  imparts  a spuriously  large  value 
to  the  attenuation  at  a line  center  at  high  altitudes  because 
the  attenuation  at  a line  center  varies  inversely  with  the  line 
width  when  Doppler  broadening  is  neglected,  while  the  pressure- 
broadened  line  width  narrows  with  increasing  altitude.  The 
effect  of  the  decrease  in  line  width  with  altitude  will  tend 
to  counterbalance  the  decreasing  concentration,  leaving  a ten- 
dency to  a spuriously  constant  attenuation  at  the  center  of  a 
line.  This  effect  has  been  allowed  for  in  the  SLAM  program  as 
described  above.  However  the  effect  of  Doppler  broadening  on  the 
wings  of  the  lines  is  not  yet  included.  This  would  involve 
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combining  the  Doppler  broadening  with  collisional  broadening  not 
only  when  the  latter  is  described  by  the  Lorentz  profile  but  also 
when  it  is  described  by  the  Van  Vleck-Weisskopf  or  Gross  profiles. 
Another  high-altitude  effect  that  must  be  considered  is  Zeeman 
splitting  of  the  oxygen  lines. 

At  lower  altitudes,  the  profile  of  two  overlapping  lines 

of  a single  molecular  constituent  is  not  necessarily  equal  to 

the  sum  of  the  profiles  of  the  individual  lines  (as  assumed  in 

SLAM) . This  consideration  is  particularly  important  for  the 

microwave  spectrum  of  oxygen.  Recent  results  suggest  how  this 

24 

can  be  taken  into  account  in  the  future. 

Finally,  we  note  the  influence  of  the  method  of  integration 
upon  the  computed  vertical  attenuations.  A rough  calculation 
in  one  particular  case  showed  an  8%  difference  between  the  re- 
sults of  the  trapezoidal-rule  integration  and  a higher-order 
method. 
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mRiTE(111a97 I VBOTj  VTOPjONU (1 l<ONU( X 1 ) * X 1 
97  RORMAT ( «VB0T'*F12«3*  * VTOP' *F12‘3* ' QNU ( 1 > ■ • * F12 *3*  '0NU-'#F12* 


3 • 


•II* 


*RlB» 

XB-1 

V2P-0NUI 11  CBOUNO 


XB-1 

P0-1013.00 

T0-29**00 


CBi(N)  AND  CB*(M,N)  ARE  FACTOR*  NECEBBARY  FOR  THE  CALCULATION  OF 
THE  POPULATION  FACTOR  ANO  THE  ROTATIONAL  PARTITION  FUNCTION  AT 
T ( n ) « RELATIVE  TO  29*  OEOREE*  KELVIN 


00  2 N-l/33 

2 CBl ( N ) ■ ( TO-T ( N ) l/(TU-T(N»-.*9A*l 
00  21  M-1,7 

00  T0(17«19«17#19«1*«17«19)M 
17  00  3 N>1#31 

3 CB2(H*N)«(  (T0/T(NII—1»BI 
00  TO  21 

19  00  A N>1*33 
A CB2 ( M/N ) «TO/T ( N I 
21  CONTINUE 


FACTOR  OIVXNO  TEMPERATURE  AND  PRkBBURE  DEPENDENCE  OR  LINE  NXOTH  IB 
COMPUTED  AT  EACH  ALUTUOE 
00  S N-l«33 

B CA(NI-( (T0/T(N))**0*9I*(P(N)/P0I 
C OOPl (N1/N2 ) II  A FACTOR  NECEIBART  TO  DETERMINE  THE  DOPPLER 

C BROAOENINO  FOR  A OIVEN  MOLECULE  AT  A OXVEN  ALTITUDE 

00  372  Nl-1#  7 
00  372  N2-l« 33 

372  OOPl (N1«N2)-00P(N1 )-T(N2l**(9 
IUM-0'0 
BAvE-0*0 
V>Vl 
C 


f 
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c 

c 


c 

c 

c 


c 

c 


M£OUENCyCOMPUTATION  °F  ATTENUAno*  AT  * PIKED  ALTITUOE  ANO 

„ OO if "'T °“r T *"’«"*•"  I"*'oETi.niNCON  'L“‘  ,0U,,° 

IF  ( V-BOuNO-flNUa  ))2»,29<  33 

29 

QO  TO  38 
33  CONTINUE 

13-11 
00  TO  83 
33  00  39  J-ISill 

37  ul!:J°UN0-,,NU'Jn37'37'3» 

00  TO  43 
39  CONTINUE 
:*-:i 

43  00  i N-l j 33 
00  *7  N-l j 7 
CAviini-o.o 
27  •uniim-o.o 
DO  23  I-ISjI* 

n-fiOL  ( 1 1 

*'  """  •'  

«-«0L  C I ) 

Z-v-qnu ( I ) 

lp-lpi 

"nlil  rSoos'cnl. I To?HIC"i ' voist  ’?I!d J!l  c°"put,°  *• 

Zl-V-ONU I I ) 


00*TO?200#201#20^Jo3^f9LATI°Ni  0EPENOI*0  ON  LINE  OHAPE 
C L0NENT2  LINE  SHAPE 

°°  ooMTof2aB**AU 

C VAN-VLECK  UEISSKOPr  LINE  ■um 

“ MKSSS^JKSSaStSS-'' 

' ~ SssSSS8*' 

00  TO  228 

c 203  »»  r«  •»**«  «.i« ,,  c.lcul.tE= 
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C THt  NEXT  PARAMETER  DEPENDS  ON  TMt  RATIO  OF  THE  COLLISIONACLT 

C SRwAOENEO  nIOTM  TO  TmE  OOPPLER  MOTH 

Yb*LPHA1/BETA 

OEPENOINO  ON  The  ABUVE  PARAMETER  EITHER  an  ASYMPTOTIC  FORMULA 
PO«  the  VOIOT  PROFILE  at  The  CEMER  OP  a RESONANCE  LINE  It  UtEO 
OR  A POLTnOMIAL  APPROXIMATION  FOR  the  PROBABILITY  INTEORAL  IB 
AOwPTED 

IFIY<QE<S< IQO  TO  20* 

TY-l./l 1<4<32759114YI 
T Yi«TY*TY 
T Y j*T Y2*TY 
TY4-TY3-TY 
TY5-TY4-TY 

8Ufll(MI-SS*l<T7245A*( <28AB2»**TY-t2BAR»*7*TY241<R2lRlR*TY3 

2-.<453152*TY4*l <0*1^0S*TY5|/BETA 
QO  TO  2£S 
204  Y2*Y«Y 
Y3bY2#Y 
Y4-Y34Y 
YS«Y44Y 
T*bY54Y 
Y7bY4*Y 
Y|bY74Y 
Y*bY|4Y 

Sunil  M I «SS4  I 1 * /T* <5/734 < 75/ YB- 1 • B79/Y74N«S4i5/Y9  I /BETA 

225  CAY1IM.I'CAY11M)4BUM1(M) 

23  COnTINc. 

CAY-0  0 
00  47  Mb  1/ 7 

47  CAY-CAY-CAV1(MI*N(H#NI 
OPU-CAY«1<3S240 
IF(N<EQ<1 IQO  TO  10 

NEl.T  THE  ATTF-  IN  aT  THE  VARIUUt  LEVELB  IB  BUMMEW 
SUn-&UM4<3*(  i N»l i l • ( CPD*BAVE I 

SO  0 i *• '"  N i IVI*Sl<  h ( Ni  I V I bOPD 

S'  .*  3P0 
4 r .NUE 
Fr,„t  IVIBV 
00  109  N-l*33 

109  UP(NjIV)bSUM-DOwN(N«IV) 

IF ( IV40VI <0T<V2PI00  TO  S3 
IFIIV4<S4DVI<0E<V2IU0  TO  S3 
IF  I I V • QE • 100  1 00  TO  S3 

IVbI V4I 

tne  frequency  is  Incremented  and  the  calculation  performeo  again 

VBVl4( IV-l I4DV 

sun-o <o 
bavE-o.o 
00  TC  28 

53  00  T0I215/214«217I«LP1 
215  *RlTE(ill*i0i>IV*V*V2P*IC 
00  TO  219 

101  FORMAT! I5*2F10*4* • LORENTZ  PROF ILE •* 3X* ' CONSTITUENT-  *181 
214  HRITE(111*22CIIV«V*V2P«IC 
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c 

c 

c 


00  TO  219 

280  FORMAT! JS/ZFIO. 4/ ' VAN  VLECK-WtXBSKOPF  PR0FILE'/3X/ 'CONSTITUENT* 
Z '/ 19 ) 

217  WRITE! 1 11*221 1 IV/V/V2P, XC 


221  FORMAT! IS/2F10*4/ ' QROBS 
112  F0RMAT(4I3/4F«*2/4X/I1/I3/ 
819  FXNALaVl+99«*0V 


PROFXLE'OX/  'CONITXTUENTa'/IBI 
' (07*2 l '/ 1QX/FS*2) 


00  107  J«1/XV 


the  next  statements  generate  an  output  file  to  be  used  as  inrut  for 
* GRAPHICS  routine  CALLEO  0RL0T1U  " 
WRITEU09/110I10/0/0/0 


110  ANO  VERTICAL  molecular  attenuation' I 

1 , , rSII.i ’ 10'0'0'0'  <ROOEL(KK/MOB)/KKal/7)/FNUI JI/BOUNO 

***  **®/***/A2#  'FREQ"'/F8*3/  * CM*4*i'/'  SOUNO"  ' / F 7 <3*  ' CM44-1 ' 


00  T0(210/211/212)/LP1 

210  WRJTEI 109/184110, 0,0,0, IC 
00  TO  21S 

1^!4»IS??,AJiiI?^0"D0NN<U,,  u"URIOS)  H»HOR  ( OB/KM  I L0RENTZ*/8Xj*  CON 
*«TXTUENT*  1 # 15 ) 

211  KRITEI109, 214110, 0,0,0, IC 

U-^'0..  H-mOR.OB/KM,  VAN  VLECK-WEXBSKOP 

00  TO  21S 

212  WRITEI109, 213110, 0,0,0, 1C 

Z*^2ITUENT« S 3/JO"OOWN(OB)  U-UPION)  M-HONIOB/KM)  QROSS'/ZX,*  CONST 

21S  WRITE(109,112)2,0,i,i,«7,7o«,(c,*2, l,3,ti 
WRITE  1 109, 180 12/90/ 1/ 1/ 1 */i000*/*2/ *2«l/8/ *1 
18C  FORMAT! 41 3/ 2FS • 0/ 8FS • 8/ 4X/ II, 13/ ’ ( OB .3  I ' / i OX/FS *2 ) 

WRITE ( 109/ 113 14/ 1/-1,0/- 1/0/0 

113  FORMAT ( 713  I 

WRITE (109/ 114  1 3/ 0/1/ *30/0 

114  FORMAT  I 81 3 I 


WRI 1 C * 109/ 1 IS  I I H ( NN ) / NN"2/ 31 1 
118  F0RMAT(10FS*0I 
1U  FORMAT (2131 

WRITE  1 109/ 181 13/90/ 1/4  30/0/ 1/1/1*/ 1000' 
WRITE! 109/122)  OOWNINN/ J)/NN«8/31 I 
122  FORMAT! 10(01,3)  | 

WRITE! 109/114 ) 8/ 42, SB/ 0 
WRITE! 109,181)3,90/1/430,0,1,1, i«,xo00* 
WRITE! 109/1881 (UR(NN/ j|/NNa2/31 I 
WRITE  1 109/ 114 1 8/42/ BB/O 


WRITE (109/ 184  1 8/ 0/1/1/ *7/70*/3/3/l/l,*l 
184  F0RMAt(4I3/4FS*2/4X/U/I3/ • ( U/-X/  ' * * ' I ' ,BX/ FS  *2  ) 
WRITE! 109/113)4, -2/0,0, 0/0,0 
WRITE! 109/ 188l2/90/l/l/*l/i00*/*Z/*2/l/7/,l 
182  FORMAT (4I3/2FS*2/2FS*2/4X/I1/I3/ • ( 4X, 07 . 2 | * , 7X/FB*8 I 
WRITE  1 109/ 113  1 4/0/ -8/ 0/0/ 0/0 
WRITE! 109/181 13/90/ 1/430,0,1,1/ *1,100* 

181  FORMAT! 713/ 8F8*2) 


WRITE ( 109/ 123  I (HORINN, J),NN*8/31  > 


183  FORMAT (10(01*3)  ) 
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WRITE!109,lUI8,-*,7*,0 
WRITE (109, 183 1 It, *,*,0,0, 931  418 
IBS  FORMAT (BI3/*FI>0,  ' Of  *ER  KM*  I 

WRITE  I 109, 117  > lt/0/1, 1/0/i 90/ IB 
117  FORMAT ( flll, 111, 'VERTICAL  HElOHT(XMI') 

WRITE<109,11|)1*,*,*,-1,0,  B/ ABB 
111  FORMAT (BIJ/tll/ 'ATTENUATION  01 ' I 
WRITE!  109/llAU 
C 

C THE  NEXT  STATEMENT!  GENERATE  A TABULAR  OUTRUT 
XRITEI 111/1 OA>FNU(JI 
10A  FORMAT! ' FREQUENCY* ' /F10 • 3 I 
WRITE ( 1 1 1/ 104  I 

104  FORMAT!*!'  HEIGHT  HOR  OOWN  UR  'll 

107  WRITE!  Ul/lOB) (H(NN)/H0R(NN/ J)/OOWNtNN/ J)/UR(NN/U)/NN»1/33I 
10B  FORMAT  1 * ( Fl 0 • 3/ 3E 1 1 • fl  I ) 

IF! ( v*»B*OV) •OE'VtiaO  TO  76 
IF! I V+DV I <QT • VtR I QO  TO  A7 
IF ( V *GE (FINAL  1 00  TO  A3 
00  TO  7B 
A3  Vl-FINAL-OV 
IB-1 

IV- 1 

V- Vl 

00  TO  IB 

NEXT  THE  OATA  FROM  THE  OATa  FILE  WILL  BE  REORBANIZBO  AND  THE  FILE 
WILL  BE  REAO  aqain 
A 7 IV-1 

VI- V-DV 
VB0T-V1-B0UN0 
00  A9  IN-1/11 

IF(ONU(IN).OT»VBOT)00  TO  71 
A9  CONTINUE 
IN-11 
71  IJ-IN 
t-1 

00  73  I-IJ/I1 
ONUlL l-ONU! 1 1 
• (L l-B! I I 
ALRMA ( L I-ALRhA ( 1 1 
EORIL l-EOR! 1 1 
MOL(LI-MOL(I» 

73  l-L*l 
I-L 
ILL-L 
00  TO  1 
7B  CALL  EXIT 
(NO 
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